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ABSTRACT 

We present a spectroscopic analysis of an extremely rapidly rotating late O-type star, VFTS102, 
observed during a spectroscopic survey of 30 Doradus. VFTS102 has a projected rotational velocity 
larger than 500 km s -1 and probably as large as 600 km s -1 ; as such it would appear to be the most 
rapidly rotating massive star currently identified. Its radial velocity differs by 40kms~ 1 from the 
mean for 30 Doradus, suggesting that it is a runaway. VFTS102 lies 12 pes from the X-ray pulsar 
PSR J0537-6910 in the tail of its X-ray diffuse emission. We suggest that these objects originated from 
a binary system with the rotational and radial velocities of VFTS102 resulting from mass transfer 
from the progenitor of PSR J0537-691 and the supernova explosion respectively. 
Subject headings: stars: early-type — stars: rotation — stars: evolution — Magellanic Clouds - 
pulsars: individual (PSR J0537-6910) 



1. INTRODUCTION 

In recent years the importance of binarity in the evo- 
lution of massive stars has been increasingly recognised. 
This aris es from most OB-ty pe stars residing in multiple 
systems (jMason et al.ll2009r ) and the significant changes 
to stellar properties that binarity can cause (see, for 
example, [P odsiadlows kTet al.l 119921 : lLanger et al.l 120081 : 
lEldridee et al.ll2011f ). 

Here we present a spectroscopic analysis of a rapidly 
rotating (?; eq sinz^ 600 km s -1 ) O-type star in the 30 
Doradus region of the Large Magellanic Cloud (LMC). 
Designated VFTS102 (jEvans et all I20TTI hereafter Pa- 
per lj3, the star is rotating m ore rapidly than any 
observed in recent large surveys (jMartavan et al.l 120061: 
iHunter et al.l l2009h and may also be a runaway. It lies 
less than one arcminute from the X-ray pulsar, PSR 
J0537-6910, which is moving away from it. 

We suggest that VFTS102 might originally have been 
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part of a binary system with the progenitor of the pulsar. 

2. OBSERVATIONS 

Spectroscopy of VFTS102 was obtained as part of 
the VLT-FLAMES Tarantula Survey, covering the 3980- 
5050A region at a spectral resolving power of 7000 to 
8500. Spectroscopy of the Ha region was also available, 
although this was not used in the quantitative analysis. 
Details of the observations and initial data reduction are 
available in Paper I. 

The spectra were normalised to selected continuum 
windows using a sigma-clipping rejection algorithm to 
exclude cosmic rays. No velocity shifts were observed 
between differen t epochs, although simulations (see, 
Sana et al. 2009) indicate that 30% of short period (less 
than 10 days) and effectively all longer term binaries 
would not have been detected. We have therefore as- 
sumed VFTS102 to be single and the sigma-clipped 
merged spectrum displays a signal-to-noise ratio of ap- 
proximately 130 and 60 for the 4000-4500 and 4500- 
5000A regions respectively. 

An 09: Vnnne spectral classification was obtained by 
smoothing and rebinning the spectrum to an effective 
resolving power of 4000 and comparing with standards 
compiled for the Tarantula Survey (Sana et al. in prepa- 
ration). The principle uncertainties arise from the ex- 
tremely large rotational broadening and significant neb- 
ular contamination of the He I lines, with the two suf- 
fixes indicating extreme line broadening ('nnn') and an 
emission-line star ('e'). 

3. ANALYSIS 

3.1. Projected rotational velocity 

The large rotational broadening of the spectral fea- 
tures makes reliable measurements of the projected ro- 
tational velocity, w cq sini, difficult. We have used 
a Fourier Transform ( FT) approach as discussed by 
iSimon-Diaz fc Herrerol ([2007D . supplemented by fitting 
rotational broadened profiles (PF) to the observed spec- 
tral features. 

The Balmer lines have significant nebular emission and 
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Figure 1. Observed spectra (in velocity space) for VFTS102 and rotationally broadened profiles for the He I line at 4026A, the combined 
profile for the He I lines at 4026, 41 43 and 4387A, the combined profile for the He II lines at 4200 and 454lA and the He II line at 4686A. 



hence the weaker helium spectra were utilized, as illus- 
trated in Fig. [Q The He I line at 447lA, although well 
observed, also showed significant nebular emission and 
was not analysed. By contrast the line at 4026A showed 
no evidence of emission and yielded a plausible minimum 
in the Fourier Transform for a w oq siniof 560 km s -1 . 
The PF methodology leads to a slightly higher estimate 
(580 km s^ 1 ). The He I lines at 4143 and 4387A were ob- 
served although they are relatively weak. They and the 
line at 4026A were converted into velocity space, merged 
and analysed. The two methodologies yielded effectively 
identical estimates of 640kms _1 ; a similar procedure was 
undertaken for the He II lines at 4200 and 454lA yield- 
ing 540 kms" 1 (FT) and 510 kms" 1 (PF). The He II line 
at 4686A was found to be sensitive to the normalisation 
with a WeqSini of ~560kms being estimated. 

The individual results should be treated with caution 
but overall they imply that this star is rotating near to 
its critical velocity, with the mean value for the FT esti- 
mates being 580 km s -1 . As discussed bv pTownsend et al.l 
(2004), projected rotational velocities may be underesti- 
mated at these large velocities. For a B0 star rotating at 
95% of the critical velocity, this underestimation will be 
approximately 10%. Hence our best estimate for the pro- 
jected rotational velocity is ^600 km s -1 . A lower limit of 
500 km s -1 has been adopted, whilst the upper value will 
be constrained by the critical velocity of approx imately 
700 km s" 1 from the models of lBrott et all (|2011l) . 



This estimate is signifi cantly higher 
(< 370 km s" 1 ) found by Martava n et al.l 



t han those 
_ (f2006h and 

IHunter et all 1)20091 ) in their LMC B-type stellar sam- 
ples. It is also larger than any of the preliminary esti- 
mates (<450 km s _1 ) for ~ 270 B-type stars in the Taran- 
tula survey, although other rapidly rotating O-type stars 
have been identified. As such it would appear to have 
the highest projected rotational velocity estimate of any 
massive star yet analysed. 

3.2. Radial velocity 

Radial velocities were measured by cross-correlating 
spectral features against a theoretical template spectrum 
taken from a grid calculated usin g the code TLUSTY 
iHubenvl (|1988l ) - see lDufton et al.l (|2005[ ) for details. 

Five spectral regions were considered, viz. H<5 and H7 
(with the cores excluded); He I at 4026A; 4630-4700A 
with strong multiplets due to C III and O II and an 
He II line; 4000-4500A (with nebular emission being ex- 
cluded). The measurements are in excellent agreement 
with a mean value of 228±12kms _1 ; if the error dis- 
tribution is normally distributed the uncertainty in this 
mean value would be 6kms _1 . 

From a study of ^180 presumably single O-type stars 
in the Tarantula survey Sana et al. (in preparation) find 
a mean velocity of 271 kms -1 with a standard deviation 
of lOkms -1 . Preliminary analysis of the B-type stars in 
the same survey has yielded 270±17kms- 1 . VFTS102 
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Table 1 

Properties of VFTS102 



Parameter Estimate 

Spectral type 09: Vnnnc 

T eff (K) ' 36000±5000 

log g (cms -2 ) 3.6±0.5 

UeqSini(km s -1 ) 600 ± 100 

■u r (kms -1 ) 228±6 

log L/L 5.0±0.2 

lies more than two standard deviations away from these 
results, implying that it might be a runaway. 

3.3. Atmospheric parameters 

While the equatorial regions of VFTS102 will have 
a lower gravity than the poles (because of centrifu- 
gal forces), and hence a lower temperature (because 
of von Zeipel gravity darkening), we first characterise 
the spectrum by comparison with those generated with 
spatially homogeneous models, convolved with a simple 
rotational-broadening function. We have use d both our 
TLUSTY grid and FASTWIND calculations ([Puis et al.l 
2005), adopting an LMC chemical composition. For the 
former, the strength of the He II spectrum implies an ef- 
fective temperature (T eS ) of -32500-35000 K, whilst the 
wings of the Balmer lines lead to a surface-gravity esti- 
mate of ~3.5 dex (cgs). For the latter after allowing for 
wind effects, the corresponding parameters are 37000 K 
and 3.7 dex. The helium spectra are consistent with a so- 
lar abundance but with the observational and theoretical 
uncertainties we cannot rule out an enhancement. 

Given its projected equatorial rotation velocity, 
VFTS102 is almost certainly viewed at smi ~ 1. Hence 
the relatively cool, low-gravity equatorial regions will 
contribute significantly to the spectrum. Although their 
surface flux is lower than for the brighter poles, the anal- 
yses discussed above may underestimate the global effec- 
tive temperature and gravity. However, the rotating-star 
models discussed below suggest that the effects are not 
very large. We therefore adopt global estimates for the 
effective temperature of 36000 K and 3.6 dex but note 
that the polar gravity could be as large as 4.0 dex. Vary- 
ing the global parameters by the error estimates listed in 
TableQ]leads to significantly poorer matches between ob- 
servation and the standard models, but, given the caveats 
discussed above, those errors should still be treated with 
caution. 

For near critical rotational velocitie s , the stellar mass 
can be estimated. iHowarth fe Smithl (|2001h show that 
the stellar mass can be written in terms of lj/lj P^. v cq 
and the polar radius. Assuming that sini ~ 1 and adopt- 
ing the critical velocities from our single star models, we 
can estimate the first two quantities. Additionally for 
any given value of uj /w c , the polar radius can be inferred 
from the absolute visual magnitude and the unreddened 
(B-V) . The former can be estimated from the luminosity 
(see Sect. 13. 4|) and the latter from our effective tem- 
perature esti mate and the L MC broad-band intensities 
calculated bv lHowarthl (|201ll ). We find M > 20 M for 
v cq ~ 600 km s" 1 and T eff < 38000 K. Only by adopting 

14 The ratio of the equatorial angular velocity to that at which 
the centrifugal acceleration equals the gravitational acceleration. 
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VFTS102. The evolutionary tracks (identified by their mass) have 
rotational velocities of approximately 600kms — 1 (dashed lines) 
and 400 kms -1 (dashed-dot line). Also shown is the evolution of 
the secondary star fol lowing mass transfer for the binary model of 
ICantiello et~aTl ll2007h (solid line). 

a smaller value for v cq can we push the mass limit down, 
but even with v cq ~ 500 km s -1 the mass must exceed 
~17M . 

3.4. Luminosity 

From extant photometry (see Paper I), the (B-V) 
colour of VFTS102 is 0.35, implying an E(B-V) of 0.6 us- 
ing colours calculated from our TLUSTY grid. Adopting 
a standard reddening law leads to a logarithmic luminos- 
ity (in solar units) of 5.0 dex, with an E(B-V) error of 
±0.1 corresponding to an uncertainty of ±0.1 dex. How- 
ever there are other possible sources of error, for example 
deviations from a standard reddening law and hence we 
have adopted a larger random error estimate of ±0.2 dex. 

As VFTS102 is an Oe-type star, its intrinsic colours 
may be redder than predicted by our TLUSTY grid 
and indeed an infrared excess is found from published 
(de-reddened) 2MASS photometry. Inspection of a K- 
band VISTA image shows no evidence of contamina- 
tion by nearby sources. Further evidence for circumstcl- 
lar material is found in the strong Ha emission, which 
is double peaked as is the nearby He I line at 6678A, 
which supports our adoption of a sini ~ 1. Addition- 
ally there are weak double-peaked Fe II emission features 
(e.g. at 4233A), consistent with an Oe-type classifica- 
tion. Unfortunately our photometry and spectroscopy 
are not contemporaneous but if VFTS102 was in a high 
state when the optical photometry was taken, we may 
have overestimated the luminosity of the central star (see 
Ide Wit et aLll2006L for colour and magnitude variations 
of Be stars). 

4. PAST AND FUTURE EVOLUTION 

Stellar evolution calculations for both single and 
binary stars are avail able in the literature (see 
iMaeder fe MevnetJ [20ll . For very fast rotation, 
they suggest that rotational mixing is so efficient 
that stars may evolve quasi-chemically homogeneously 
(IMaeder I [19871: IWooslev fe He eed [20061: ICantiello et all 
120071: Ide Mink et all 120091: iBrott et al.l 1201 If ). However, 
with different physical assumptions, models do not evolve 
chemically homogeneously even for the fastest rotation 
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rates (|Cantiello et al.l 120071; lEkstrom et al.ll2008D . 

4.1. Single star evolution 

Fig. [5] illustrates evolutionary tracks for LMC single 
stars calculated using the methodology of iBrott et aD 
(|2011[) for an initial equatorial rotational velocity of 
600 km s -1 , together with that for a more slowly rotat- 
ing model. The former are evolving chemically homo- 
geneously wdiUs^rtifiJat^erJoJlows a 'normal' evolution- 
ary path. lEkstrom et al.l (|2008f) calculated models for a 
range of metallicities and masses between 3 and 60 M© 
but found that the stars followed normal evolutionary 
paths even for near critical rotational velocities. 

The estimated parameters of VFTS102 are consistent 
with our tracks for initial masses of ^20-30 M Q . Our 
models show a relatively rapid increase in the surface he- 
lium abundance due to their homogeneous evolution. For 
example the 25 M model shows an enrichment of a fac- 
tor of two after approximately 4 million years and when 
the effective temperature has incre ased to approximately 
39000 K. By contrast the models of lEkstrom et al.l (|2008f ) 
show no significant helium abundance implying that an 
accurate helium abundance estimate for VFTS102 would 
help constrain the physical assumptions. 

4.2. Binary star evolution 

Below, we first discuss the environment of VFTS102 
and then consider a possible evolutionary scenario. 

4.2.1. A pulsar near VFTS102 

VFTS102 lies in a complex environment near the 
open cluster NGC 2060. In particular it lies close to 
a you ng X-ray pulsar PSR J0537-6910 ()Marshall et al.l 
1998) and the Crab-li ke supernova remnant B0538-691 
(jMicelotta et al J 120091 ) . VFTS102 has an angular sepa- 
ration of approximately 0.8 arcminutes from PSR J0537- 
6910 implying a spatial separation (in the plane of the 
sky) of approximately 12 pc. 

The X-ray emission consists of a pulsed localised com- 
ponent and a more spatially diffuse component, with the 
latter providing the majority of the energy. The dif- 
fuse component was identified in ROSAT and ASCA 
observations by iWang fe Gotthelj (|1998a|) and inter- 
preted as coming from ram-pressure-confined material 
with the X-ray pulsar b eing identified soon afterwa rds by 
IMarshall et"aU (fl998l) . IWang fc Gotthelfl (|1998bl) anal- 
ysed ROSAT HRI observations and suggested that the 
emission could come from the remnants of a bow shock 
if the pulsar was m oving with a velocity of ^lOOOkms -1 . 
IWang et al.l (|2001| ) subsequently analysed higher spatial 
resolution CHANDRA observations, which clearly de- 
lineated this emission and implied that the pulsar was 
moving away from VFTS102. See, however, Chen et 
al. (2006) for an alternative explanation of the X-ray 
morphology, which does not imply a high velocity for 
the pulsar. Fig. [3] superimposes these emission contours 
onto an HST optical image with VFTS102 being near 
the ta il of these contours. As discussed by IWang et al.l 
(|2001D the spatial distribution of the diffuse X-ray emis- 
sion and the SNR optical emission are well correlated. 
Differences probably arise from a foreground dark cloud 
and photoionization and mechanical energy input from 
the nearby open cluster. 
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Figure 3. HST WFPC2 V-band (F606W filter) image with con- 
tours from the smoothed Chandra HRC-I image overlaid. The 
positions of VFTS102 and PSR J0537-6910 are labeled. 

Timing measurements im ply that the pulsar ha s a char- 
acteristic age of 5000 years (IMarshall et al.lll998f) , consis- 
tent with the age estimate of IWang fc Gotthelfl (IT 998b) 
from analysis of X-ray emission. iSpvrou fc Stergioulasl 
(2002) discuss the estimation of ages from spin rates and 
find the results to be sensitive to both the breaking index 
and the composition of the pulsar core. Indeed phase 
connecte d braking index m easurements for young pul- 
sars (see lZhang et al.l[200lL and references therein) yield 
breaking indices lower than the n=3 normally adopted 
with corresp onding increases in the characteristic ages. 
Additionally, iChu et al.l (|1992l) found an age of approxi- 
mately 24000 years from the kinematics of the supernova 
remnant. 

Adopting an age of 5000 years would imply that if 
these objects had been part of a binary system, their 
relative velocity (v s ) in the plane of the sky would be 
approximately 2500kms _1 . Increasing this age to 24000 
years would then imply v s ~ 500kms _1 . These values 
although large are consistent with a pulsar ve l ocity o f 
1000 kms- 1 in the model of IWang fe Gotthelfl (Il998bl ) 
and of ~600kms _1 from the separation of the diffuse 
X-r ay and radio e missio n (jWang et al.ll2001| ). Addition- 
ally Hobbsiiil] (120051) found a mean space velocity of 
approximately 400 km s -1 for a sample of young pulsars 
with velocities as hi gh as 1600 km s" 1 . From the theo- 
retical point of view, IStonl (fl982h found supernova kick 
velocitie s normally in excess of 300kms _1 , while more 
recently lEldridge et al.l (|2011f ) estimated kick-velocities 
for a single neutron star of more than 1000kms _1 with a 
mean value of ^500 km s -1 . 

4.2.2. A binary evolution scenario for VFTS102 

While the fast rotation of VFTS102 might be the result 
of the star formation process, it could also have arisen 
from spin-up due to mass transfer in a binary system 
(|Packetlll981| ). A subsequent supernova explosion of the 
donor star could the n lead to an anomalous rad ial ve- 
locity for VFTS102 (|Blaauwl fl96lT IStonel fl982l) . The 
nearby pulsar and supernova remnant make this an at- 
tractive scenario. Of course, we cannot eliminate other 
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possi ble scenarios, e.g. dynamical eject ion from a cluster 
(see iGvaramadze fc Gualandrisll2011l) but it is unclear 
whether these could produce the very large rotational 
v elocity of VFTS 1 02. 

ICantiello et al.l (|2007j ) have modelled a binary system 
with initial masses of 15 and 16 M Q adopting SMC metal- 
licity. After mass transfer the primary exploded as a type 
Ib/c supernova. At that stage the secondary has a mass 
of approximately 21 M , a rotational velocity close to 
critical and a logarithmic luminosity of approximately 
4.9 dex (see Fig. [2] for its subsequent evolution). These 
properties closely match the estimates for VFTS102 sum- 
marized in Table [T] 

Based on grids of detailed binary evolutio nary models 
(jWellstein et all 120011: Ide Mink et all l2007f ) . the initial 
masses of the two components of such a binary system 
should be comparable, with M2/M1 > 0.7. If the ini- 
tial mass of the secondary was in the range of 14-18 M Q , 
that of the primary would need to be smaller than about 
25 M Q . This agrees with the estimated initial mass of 
the supernova proge nitor based on the kin ematics of the 
supernova remnant (jMicelotta et al.ll2009| ). 

In this scenario, it takes the primary star about 11 Myr 
to evolve to the supernova stage. While the most mas- 
sive stars in 30 Dorad us have ages of a few million years 
ijWalborn et al.l fl~999f ). there is also evidence for differ- 
ent massive stella r populations with ages ranging up 
to about 10 My r (IWalborn fc Blades! I1997D . Recently, 
iDe Marchi et al.l ([20111 ) have undertaken an extensive 
study of lower mass (^4 M Q ) main sequence and pre- 
main sequence stars in 30 Doradus. They obtain a me- 
dian age of 12 Myr with ages of up 30 Myr. Hence 
it would appear possible that the putative binary sys- 
tem formed in the vicinity of 30 Doradus approximately 
10 Myr ago and underw ent an evolutionary hi story sim- 
ilar to that modelled bv ICantiello et all (|2007l ). 

While the evolutionary link between VFTS102 and 
the pulsar is an attractive scenario, the offset between 
VFTS102 and the apparent center of the radio emis- 
sion (see Wang et al. 2001) remains to be explained. 
Proper motion information would be extremely valuable 
to further test this hypothesis. PSR J0537-6910 has not 
been definitely iden tified in other wavelength regions. 
iMignani et all (2005) using ACS imaging from the Hub- 
ble Space telescope found two plausible identifications 
that would imply an optical luminos ity similar to the 
Crab- like pulsars. A radio survey by iManchester et all 
(pOOel only yielded an upper limit to its luminosity con- 
sistent with other millisecond pulsars. However esti- 
mates for both components may be obtained from the 
HST proper motion study (Programme: 12499; PI: D.J. 
Lennon) that is currently underway. 

4.3. Evolutionary future 

Irrespective of the origin of VFTS102, it is interesting 
to consider its likely fate. Stellar evolutionary models 
of rapidly rotating stars h ave recently been gen erated 
bv lWooslev fc Hegerl (pOOfih and lYoon et all pM). The 
latter consider the fate of objects with rotational veloci- 
ties up to the critical value (v c ) . The evolution is shown 
to depend not only on initial mass and rotational veloc- 
ity but also on the metallicity. In particular GRBs are 
predicted to occur only at sub-solar metallicities. 

Based on our single star models, VFTS102 has a 



rotational velocity above ~ 0.8v c and is thus ex- 
pected to ev olve quasi- chemically homog eneously. While 
lYoon et al.l ((20061) and lWooslev fc Hegerl (120061) estimate 
the metallicity threshold for GRB formation from chem- 
ically homogeneous evolution to be somewhat below the 
LMC metallicit y, the latter note its s ensitivity to the 
mass loss rate (|Vink fc de Koterll2005D . Indeed all our 
most rapidly rotating 20 — 30 M models are evolv- 
ing chemically homogeneously throughout core hydrogen 
burning (Fig. [2]), a prerequisite to qualify for a GRB 
progenitor. In any case, within the context of homoge- 
neous evolution VFTS102 is expected to form a rapidly 
rotating black hole, and a Type Ic hypernova. This con- 
jecture remains the s ame within the binary scenario of 
ICantiello et all ((2007I) . 

Assuming a space velocity of 40kms~ 1 for VFTS102 
(compatible with its anomalous radial velocity), our 
evolutionary models imply that VFTS102 will travel 
~300-400 pc before en ding its life. This is consis- 
tent with the finding of iHammer et all (j2006l ) that the 
locations of three nearby GRBs were found several 
hundred parsecs away from their most likely progeni- 
tor birth locations (see, however, iMargutti et al.l 120071 : 
iWiersema et al.ll2007t lHan et al.ll20Tfj[l 

5. CONCLUSIONS 

VFTS102 has a projected rotational velocity far higher 
than those found in previous surveys of massive stars in 
the LMC, and indeed it would appear to qualify as the 
most rapidly rotating massive star yet identified. With 
a luminosity of 10 5 L Q we estimate its current mass to 
be approximately 25 M Q . Its extreme rotation, peculiar 
radial velocity, proximity to the X-ray pulsar PSR J0537- 
6910 and to a supernova remnant suggest that the star 
is the result of binary interaction. 

It is proposed that VFTS102 and the pulsar origi- 
nated in a binary system with mass transfer spinning-up 
VFTS102 and the supernova explosion imparting radial 
velocity kicks to both components. If evolving chemically 
homogeneously, as suggested by recent models, VFTS102 
could become a GRB or hypernova at the end of its life. 
Additionally it may provide a critical test case for chem- 
ically homogeneous evolution. 
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